Drift-Alfven vortex structures in the edge region of a fusion relevant plasma 
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Edge turbulent structures are commonly observed in fusion devices and are generally believed 
to be responsible for confinement degradation. Among their origin Drift-Alfven turbulence is one 
of the most commonly suggested. Drift-Alfven paradigm allows the existence of localized vortex- 
like structures observed also in various systems. Here we present the evidence of the presence of 
drift- Alfven vortices in the edge region of RFX-Mod RFP device, showing how these structures are 
responsible for electromagnetic turbulence at the edge and its intermittent nature. 



Turbulence represents an outstanding critical issue in 
the physics of magnetically confined plasmas for ther- 
monuclear fusion research. Indeed plasma turbulence has 
been recognized since the beginning as the cause of the 
so-called anomalous particle and energy transport 
In recent years it has been observed that, within inco- 
herent fluctuations, coherent structures emerge similar 
to vortices observed in fluid turbulence These struc- 
tures have been detected in a variety of devices, ranging 
from tokamaks 0, 0, 0, 0L through stellarators 0] , up 
to reversed field pinches [81] and linear devices [9(], and 
represent a features shared with astrophysical plasmas 
|ld |. Turbulent structures, often referred to as blobs, are 
responsible for the high degree of intermittency gener- 
ally observed. Indeed the generation of these structures, 
arising because of the presence of various instabilities in 
non-linear regime, is responsible for the breaking of self- 
similarity in the energy cascade process [11 1 



Blobs arising in fusion relevant plasmas, have been ex- 
tensively studied in the plane perpendicular to the main 
magnetic field and only recently their 3D features 
have been experimentally addressed [I4I . This interest is 
enhanced by some analogies with Edge Localized Modes 
(ELMs), which arc indeed thought to be associated with 
parallel current filaments [l3| ■ Present theories about 
blob formation and dynamics suggest an interchange-like 
origin, with effects induced by sheath boundary condi- 
tions of the material objects intersecting the magnetic 
flux surface [l^ . Plasma quasi- neutrality implies the con- 
dition V • J = on the total current J: considering the 
non vanishing V_l components of the diamagnetic and 
polarization curreiit, a parallel current density perturba- 
tion j'li must arise [15| . Experimental evidence on the ex- 
istence of filaments associated with blobs have been found 
Nevertheless, although interchange is believed to 
be responsible for blobs in the Scrape Off Layer plasmas, 
they do not represent the only possible mechanism for the 
generation of electromagnetic coherent structures. Drift 
wave instability, which is thought to dominate plasma 
turbulence in the edge region |l6| . is a nonlinear, non 
periodic motion involving disturbances on a background 
pressure gradient of a magnetized plasma and eddies of 
fluid-like motion in which the advecting velocity of all 
charged species is the E x B velocity [17[. On the theo- 



retical level the electron dynamics is purely electrostatic 
if the parallel Alfven transit frequency is much faster than 
the thermal electron transit frequency (or equivalently if 
f3 me /Mi) and large than any drift- wave frequency. 
Actually these conditions are not satisfied in the edge 
region of fusion devices [l6| and the resulting turbulence 
and transport level will be determined by electromagnetic 
effects in the framework of drift- Alfven dynamics, which 
represents the paradigm for the description of the cou- 
pling of drift-waves with Kinetic Alfven waves (KAW) 
[l8|. The key distinction between drift- Alfven dynam- 
ics and magnetohydrodynamics (MHD) is the inclusion 
of parallel electron motion and electron pressure effects. 
The disturbances in the electric field arising from the 
presence of fluid eddies are caused by the tendency of 
the electrons to establish a force balance along the mag- 
netic field lines. Pressure disturbances have their paral- 
lel gradients balanced by a parallel electric field, whose 
static part is given by a parallel gradient of the electro- 
static potential. Turbulence itself is driven by the back- 
ground gradient and the electron pressure and electro- 
static potential are coupled together through parallel cur- 
rents. The corresponding magnetic fluctuations could not 
contribute to a direct enhancing of cross-fleld transport 
through the so-called magnetic flutter transport [31 j but 
provide an additional coupling between parallel drift cur- 
rent and electrostatic drift wave potential As shown 
in electromagnetic effects are important for drift- 
wave dynamics at much smaller values of plasma /3 than 
expected by pure MHD considerations. This causes drift- 
wave dynamics to compete with the effects of the rather 
small parallel electron resistivity, the latter responsible 
(in the drift-wave framework) for the phase relationship 
between density and potential fluctuations [Tij . 

In the non- linear regimes, drift- Alfven turbulence may 
generate non linear structures in the form of electromag- 
netic vortices [l^l- As aforementioned these structures 
are generated by the non linear coupling of drift waves 
and Kinetic Alfven waves. The latter exhibit a dispersion 
relation modiflcd with respect to the shear Alfven waves, 
with a term including the perpendicular wave vector , 
= [i_(_(fc^^p^)2]i/2 where va is the Alfven velocity and 
Ps the ion sound gyroradius ps = Cs/^i- These struc- 
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Figure 1: Average coherent structure detected at scale r = 
4/xs using the electron pressure as reference signal. All wave- 
forms represent variations with respect to average values: (a) 
electron pressure, (b) electron density, (c) electron tempera- 
ture, (d) plasma potential, (e) parallel current density 



tures have been observed both in astrophysical plasmas 
[13, EH and in linear devices [9[, but up to now they 
have not been experimentally observed in fusion relevant 
devices. In this Letter we present a clear experimental 
evidence of the existence of Drift- Alfvcn-Kinetic (DKA) 
vortices in the edge region of the RFX-mod reversed field 
pinch experiment 22 1. 

Despite the peculiar magnetic topology, the edge re- 
gion of RFP plasmas shares many features with other 
magnetic devices, among which the strong intermittent 
character of electrostatic fluctuations and the highly 
sheared E x B flow detected at the edge [l^. Inter- 
mittency manifests itself as a clear departure from self- 
similarity, and can be imputed to the presence of orga- 
nized structures, intermittent structures, which make the 
process of energy cascade inhomogeneous. They have 
been extensively studied from the experimental point of 
view, observing their vortex-like shape on the perpen- 
dicular plane with an associated pressure perturbation 
[1, Ell • They contribute to the cross-field transport for 
up to 50 % of the particle losses [1] . Recently their elec- 
tromagnetic features have been experimentally described, 
revealing the existence of a parallel current density fluc- 
tuation j'li associated to the pressure perturbation, which 
can represent up to few % of the total parallel current 
[l^ , but no theoretical explanation has been proposed to 
determine the underlying instability mechanism respon- 
sible for the formation of these structures. Although re- 
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Figure 2: Top: electron pressure as a function of time and 
radial position normalized to ion sound gyroradius ps- Bot- 
tom: toroidal magnetic field as a function of time and toroidal 
coordinate normalized to ps- 



sistive interchange are expected to be unstable in RFPs 
24| . Drift- Alfven turbulence could also play a role, con- 
sidering the pressure gradient and the P condition en- 
countered. 



The data presented hereafter have been obtained in the 
RFX-mod reversed field pinch device {R/a = 2m/0.459TO 
) 1 operating at relatively low plasma current (Ip < 
400kA) and with average density normalized to the 
Greenwald density n/ug « 0.4 — 0.5. The typical plasma 
parameters observed in the edge region for this type of 
discharge are density of the order of 1 — 2 x 10^^ 
m~'^, temperature in the range 20-40 eV and magnetic 
field Bq around 0.15 T. The corresponding (3 are in the 
range of 1-2 % thus ensuring the condition f3 ^ me/ Mi 
whereas the typical scale length pa = Cs/Qi is equal to 3- 
4 mm. It is worth to remember that in the edge region of 
RFP plasmas the magnetic field is essentially poloidal, so 
that the perpendicular plane corresponds to the radial- 
toroidal plane. 

A new insertable probe, developed in order to s tudy 
electromagnetic turbulence and described elsewhere [2^ , 
has been used to explore the last 5 cm of the plasma col- 
umn. The system consists of two boron nitride cases, 
each of them housing 5x8 electrostatic pins radially 
spaced by 6 mm. The pins are used as 5 pins balanced 
triple probe, allowing the simultaneous measurement of 
plasma density, electron temperature, electron pressure, 
plasma potential and their radial profiles at the same 
toroidal location, as well as the radial and toroidal com- 
ponents of the E X B plasma velocity. The particular ar- 
rangement of electrostatic measurements allows a direct 
estimate of the local fluctuation of vorticity w = V x v, 
where v is the electric drift velocity, from the floating po- 
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Figure 3: (a) Hodogram of the E x B plasma velocity in the 
perpendicular plane, showing the closed path corresponding 
to a parallel vorticity perturbation, (b) Hodogram of the mag- 
netic field fluctuations in the perpendicular plane, showing a 
closed path corresponding to the transit of a current density 
filament. 



Figure 4: Velocity and magnetic field profiles for the average 
coherent structure detected at t — 4/xs. (a) Toroidal compo- 
nent of the E X B plasma velocity, and of the corresponding 
component of Alfven velocity as computed from the magnetic 
fluctuation measurements, (b) The same of panel (a) but for 
the radial component. 



tential ones Vf, as ojjj = ^V^^V/, where plasma potential 
has been approximated by floating potential as usually 
done [i^. A radial array of 7 three-axial magnetic coils 
is located in each case, in order to measure the fluctua- 
tions of the three components of the magnetic field. Thus 
a direct estimate of the parallel current density can be 
done from Ampere's law ~ jg — -^{d^br — drb^), vir- 
tually in the same toroidal position of the vorticity mea- 
surements. Data were digitally sampled at 5 MHz with 
a minimum bandwidth of 700 kHz. The data collected 
with the insertable probe have been completed with mea- 
surements obtained from a toroidal distributed array of 
magnetic pick-up coils located inside the vacuum vessel, 
pertaining to the ISIS system 27|. The data analysis 
technique used to disentangle coherent structures from 
the turbulent background is based on wavelet analysis 
and has been extensively described elsewhere [1^. It 
allows to locate within the signal the presence of struc- 
tures at a given temporal scale. This method has been 
used together with the traditional conditional averaging 
technique to better extract the common features of the 
observed structures. 

In figure [1] the results of a conditional average pro- 
cedure are shown. All the time windows used for the 
average have been chosen using the appearance of an in- 
termittent structure at characteristic time scale t = 4/is 
(well above the l/fi,; time scale) on the pressure signal. 
It can be easily recognized in panel (a) the pressure peak, 
typical of plasma blobs, mainly determined by electron 
density (b). The resulting temperature structure display 



a doubly-peaked pattern, whose impact on the plasma 
pressure is less important, but contributes in determin- 
ing the plasma potential pattern, shown in panel (d). It 
can also be easily recognized that electron density and 
plasma potential are nearly in phase, with a phase differ- 
ence around 1 /US, suggesting the drift-origin of the ob- 
served non-linear structure. The typical pattern of cur- 
rent density associated to these measurements is shown 
in panel (e) of the same figure. It can be easily recog- 
nized the existence of a current peak associated to the 
electron pressure blob. The slight phase shift observed 
may be imputed to the small deviation of the nominal 
toroidal position of current measurement with respect 
to the pressure one, and is consistent with the E x B 
toroidal propagation of the structure [25| . The radial ar- 
ray of electron pressure measurements allows to investi- 
gate the radial extension of this pressure perturbation, as 
shown in the upper panel of figure [2] The pressure peak 
exhibits a radial extent of 2-3 ps, which is indeed the 
typical extent of the DKA vortex predicted for example 
in [iO] and observed in [13]. In the other perpendicu- 
lar direction the dimension may be estimated using the 
toroidal distributed array of pick-up coils, and examin- 
ing the magnetic footprint of the structure. The results 
of a conditional average procedure still with the appear- 
ance of an intermittent structure on electron pressures, 
is shown in the lower panel of figure [21 In the direction 
of the flow the structure is larger (of the order of 100 
Ps corresponding to approximately 40 cm). This larger 
dimension may be imputed to the stretching effect in- 
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Figure 5: Comparison between parallel current density (solid 
curve) and parallel vorticity (dashed curve) for average coher- 
ent structure detected at scale r = 4/is. 



previous figure. The two quantities are found to be very 
well correlated one to the other. This last observation 
together with the drift-type phase relation between po- 
tential and density, and the alfvcnicity of the fluctuating 
velocity, establishes without ambiguity the Drift-Kinetic 
nature of the intermittent structures observed at the edge 
of Reversed Field Pinches. These measurements suggest 
the necessity to complete blob description with a full elec- 
tromagnetic characterization, and support the theory of 
Drift- Alfven dynamics as a paradigm for the description 
of the edge confined region of thermonuclear plasmas. 
This work was supported by the Euratom Communi- 
ties under the contract of Association between EU- 
RATOM/ENEA. The views and opinions expressed 
herein do not necessarily reflect those of the European 
Commission. 



duccd by the sheared plasma flow, which is well known 
to modify vortex structures [2^. It has been already re- 
ported elsewhere that these structures are convected by 
the mean E x B flow and travel in the electron diamag- 
netic direction [25j . Using a poloidal array of pick-up 
coils pertaining to the same ISIS system it is possible to 
infer the poloidal wave vector kg w fcy of these structures 
which is of the order of 3.3 m~^. 

To conflrm the vortex nature of the observed struc- 
ture, the corresponding average vexb velocity fluctua- 
tions have been calculated. Results arc shown in figure [3] 
panel (a) where the hodogram of the drift velocity com- 
ponent on the perpendicular plane is shown. The fluctu- 
ating velocities follow a closed trajectory on the perpen- 
dicular plane, as a consequence of the vortcx-likc nature 
of the observed structure. Similarly the hodogram of the 
perpendicular components of the magnetic field (figure 
[3] panel (b)) shows a similar closed path, corresponding 
to the transit of the current density filament observed in 
figure [T](e). In figure |4] panel (a) and (b) the two com- 
ponents of the E X B velocity fiuctuations are compared 
to the corresponding components of Alfven velocity fluc- 



measurements, 



and 



The good 



match observed highlights the Alfvenic nature of the fluc- 
tuating velocities, reinforcing the hypothesis of DKA as 
underlying physical mechanism. As a final confirmation 
of the nature of the observed structures we perform a 
direct comparison between vorticity and parallel current 
density. In the drift-Alfven framework electrostatic po- 
tential Vp and parallel component of the vector potential 
A|| are intrinsically related and almost proportional one 
to the other [l^j H^- As aforementioned the potential 
measurement arrangement in the U-Probe allows the di- 
rect estimate of the vorticity as W|| = "g^l^/ and this 
can be compared to the parallel component of the current 
density jy = V^Ay. This comparison is shown in figure 
[SI the patterns of j|| and lo result from the conditional 
averaging procedure with the same condition used for the 
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